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ABSTRACT __ 

Recent observations of the Crab Nebula ( Rudy et al.||2015 ) have maintained its repu¬ 
tation for high energy astrophysical enlightenment and its use as a testbed for theories 
of the behaviour of magnetized, relativistic plasma. In particular, new observations of 
the inner knot located 0.65" SE from the pulsar confirm that it is compact, elongated 
transversely to the symmetry axis and curved concave t owards the pulsar. 60 percent 
polarization has been measured along the symmetry axis (Moran et al. 2013). The knot 
does not appear to be involved in the gamma ray flares, 'i'he new observations both 
reinforce the interpretation of the knot as dissipation of the pulsar wind at a strong 
shock and challenge the details of existing models of this process. In particular, it is 
argued that the compactness, high polarization and curvature are difficult to reconcile 
with simple relativistic shock models. Alternative possibilities include deflection of the 
outflow ahead of the shock and spatial variation in which the knot is interpreted as a 
caustic. Some future observations are proposed and new theoretical investigations are 
suggested. 
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1 INTRODUCTION 

The Crab Nebula has a prominent role in the history of 
high energy astrophysics. Ever since its discovery as a su¬ 
pernova 961 years ago, it has exhibited many of the earliest 
and the best examples of important high energy astrophys¬ 
ical phenomena including synchrotron radiation, spinning, 
magnetized neutron stars, relativistic outflows and, most re¬ 
cently, dramatic gamma ray flares (e.g., |Hester|2008| |Buhler| 
|fc Blandford|[^14 1. 

However, the details of this general interpretation re¬ 
main controversial. There should be a conversion of elec¬ 
tromagnetic to particle energy flux in the outflow and the 
extent to which this happens, as measured by the magneti¬ 
zation parameter a = Poynting flux/particle energy flux, is 
quite uncertain (e.g., Kirk et al.|2009| Arons|2012 1. Another 
uncertainty is the variation of the wind power with pulsar 
colatitude, 9. This is generally supposed to decrease from 
equator to pole but the actual form may depend on the in¬ 


clination angle of the pulsar magnetic axis (Tchekhovskoy 
et al.|20T5 I. A third issue is the fate of the magnetic stripes 


expec ted to be present in the equatorial outflow (Michel 
19711. These have a wavelength of only 10^ km, but there 
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may be insufficient time in the outflowing reference frame 
to dissipate them (e.g., |CoronitI||1990[ [Lyubarsky fc Kirk| 


20011. A final uncertainty is the nature of the deceleration 
of the outflow from ~ c to ~ 0.005c, the nebular expansion 
speed. It has been commonly presumed that this occurs at 
a single strong shock with post shock pressure balanced by 


some ambient nebular pressure (Rees & Gunn 1974 Ken- 
jnel fc Coroniti|1984p . However, if a is high and/or the first 
shock is oblique, the post-shock flow will be high speed or 
even supersonic and additional shocks may be present. The 
discovery of the inner knot ( Hester et al.|[T9^ l was there¬ 
fore of extreme importance as it promised to address these 
questions. 

The knot appears as a compact emitting feature lo¬ 
cated 0.65" southeast of the pulsar. It lies on the rotation 
axis and is elongated perpendicular to the axis, with a size 
0.32" X 0.16" in optical and 0.40" x 0.32" in infrared (IR), 


and is found to be concave towards the pulsar ( Rudy et al. 
|2015p . The emission from the knot is highly polarized, with a 
polarization degree ~ 60% and position angle approximately 


along the symmetry axis ( Moran et al.|2013 1 . This is consis¬ 
tent with a coherent toroidal magnetic field at the emission 
site. Meanwhile, the knot emission shows strong variability 
on a time scale of months to years, though no convincing 
association with the gamma ray flares ( [Rudy et al.||2015| . 

The knot has been convincingly interpreted as Doppler- 
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2 Y. Yuan & R. Blandford 


boosted emission from behind a termination shock (Komis¬ 


sarov & Lyubarsky 2004 Komissarov & Lyutikov 20111. 


However, these new observational results on the compact¬ 
ness, shape and polarization of the knot challenge this in¬ 
terpretation. 

In the following we examine the shock model carefully 
to exhibit these problems, and suggest possible solutions to 
them. In particular, summarizes the problem with the 
knot size and its offset from the pulsar, illustrates the 
curvature problem, ^discusses the Doppler-depolarization 
problem, ^presents a detailed model of the knot intensity 
profile and polarization, suggesting additional ingredients 
needed to fix these problems. In !|^ and 0 we discuss the 
spectral and temporal features of these observations which 
provide fresh clues as to the nature of the inner knot. We 
summarize our conclusions, and suggest some new possible 
interpretations of the inner knot which point to additional 
observational and theoretical investigations in ^ 

After most of the research reported here had been com¬ 
pleted, the authors learned of a complementary investigation 
by Komissarov, Lyutikov and Forth which will be reported 
elsewhere. 


2 COMPACTNESS PROBLEM 


Consider a simple adiabatic, ideal MHD shock ( 

Komissarov 

j&c Lyubarsky|2004||Komissarov & Lyutikov|2011 

1 (FigurejH. 


The jump conditions for such a shock assuming an isotropic 
pressure can be simply obtained by Doppler boosting the 
jump conditions for a perpendicular shock ( [Kennel fc Coro-] 
niti |1984 1 into a frame moving along the shock, or equiva¬ 


lently, following Komissarov & Lyutikov (2011 1 who showed 


that when the upstream flow is cold and extremely relativis¬ 
tic, the compression ratio across the shock only depends on 
the upstream magnetization parameter a = Bf fA-KnijlmR, 
where Bi, ni and 71 are the upstream magnetic field, parti¬ 
cle proper density and wind Lorentz factor, respectively (we 
use subscript 1 for upstream variables and 2 for downstream 
ones): 


X = 


V 2 n _ B\ _ I + 2a + -I- 16a -|- 1 

Vi„ B 2 6{l + a) 


( 1 ) 


where Ui„ and V 2 n are upstream and downstream flow ve¬ 
locity components parallel to the shock normal, respectively. 
It is easily seen that when a —> 0, y —> 1/3-1-4(j/3; when 
a —>■ 00 , X 1 “ l/(2o'). The downstream bulk Lorentz 
factor is 

CSC 5i 


72 




( 2 ) 


where (5i is the angle between the upstream velocity and the 
shock surface. If the shock is highly oblique and/or the mag¬ 
netization is quite high, the downstream flow would remain 
relativistic. The emission would then be highly beamed, and 
the knot corresponds to the spot where the outflow velocity 
is directed toward us. 

We note that as the shock only reduces the perpen¬ 
dicular velocity, the flow will be deflected from the radial 
direction. If we denote the angle between the shock surface 
and downstream velocity as S 2 , it can be shown that 



Figure 2. Contour plot of deflection angle A on the a —Si plane. 
The deflection angle 

A = — (52 = (5i — arctan(xtan5i) (4) 

is thus determined by a and 5i. The geometry in Figure 
indicates that the measured knot-pulsar angular separation 
should be d'p = r^A/D, where is the actual distance be¬ 
tween the pulsar and the emission site at the shock surface, 
and D ~ 2 kpc is the distance of the Crab Nebula. Typically 
= 0.65", and assuming Vk ^ req where Veq = 11.3"D is 


the radius of the shock in the equatorial plane (Weisskopf 
et al.|| 2012 | , we get A ^ 0.057 rad. From the contour plot 


of A on the a — 5i plane (Figure]^, we see that cr must be 
less than 4 for the deflection to be this large. 

On the other hand, the knot size perpendicular to the 
symmetry axis, characterized by its full width at half max¬ 
imum r]±, can be determined from the outflow Lorentz fac¬ 
tor, taking into account azimuthal symmetry. Based on 
a rough estimation that the Doppler beaming angle is 
~ 1 / 72 , we have ri± « 2rkj{'y2D). The ratio J^u/Ikp = 
2/(72 A) therefore only depends on a and and constrains 
these two parameters. We find that rj±/^p = 2 /( 72 A) = 
2 sin — arctan(xtan5i)] ^ 2 . 8 , with the min¬ 

imum happening at (j = 0 and = 0.39. However, ri± ~ 
0.54'p is observed (Rudy et al. 20151 and so there’s clear 


tan52 = xtan(5i. 


( 3 ) 


inconsistency. 

The problem remains with more complex shock models, 
for example, if the upstream plasma is not cold, then as long 
as 71 ^ 1, we obtain the same compression ratio as Equation 
([^ with the definition of a changed into a — Bif4TTWiji, 
where Wi = pi + Pi is the upstream enthalpy density. This 
does not change 72 or A. 

Alternatively, if the plasma is anisotropic though still 
gyrotropic, namely, the particle velocity perpendicular to 
the local magnetic field is isotropic but differ from parallel 
velocity distribution, then the compression ratio is given by 

2 + 3(7 -(- (/(j -(- {2 -\- 5(7 + 3(((7)^ -(- 8(j(l -\- 


X = 


4(2+ 0(1+ < 7 ) 


( 5 ) 
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Implications of Crab Inner Knot 3 


Plane of the sky 



Figure 1. Geometry of the termination shock on the meridional plane determined by the pulsar rotation axis and our line of sight toward 
the pulsar. The pulsar is located at point P. By convention, the pulsar rotation axis, when projected on the plane of the sky (the plane 
perpendicular to the line of sight, and therefore perpendicular to the plane of the paper), has a position angle = 300° north through 
east. The angle between line of sight and pulsar rotation axis is 0ob = 27r/3 (e.g., |Ng fc Romani|2004| [Weisskopf et al.|2012[ | . The knot 
is projected on the plane of the sky at an angle ikp Ri 0.65” southeast of the pulsar (represented by point K here). Other symbols are 
defined in the text. 


where cr = BljAnwi'yl = + Pix) 7 i, Px and 

P|| are the pressure tensor components perpendicular and 
parallel to the local magnetic field, respectively, and we 
assumed P 211 = C,P 2 _l (for details see Appendix All. In 
the limit ^ = 0 (no parallel pressure), we get x = (2 + 
3(t + V4 + 28cr + 25 (t^)/(8 + 8 ct). x 1/2 + 3(t/ 4 when 
(7 —> 0 and X 1 ~ 2 /( 5 ( 7 ) when a —> 00 . The ex¬ 
pressions for 72 and A stay the same with x changed to 
this new value, so we get the ratio jjx/vl'p = 2 /( 72 A) = 
2 sin (5i — X^/[<5i — arctan(xtan5i)] 3.46, the mini¬ 

mum occurring when <7 = 0 and (5i = 0: this is worse 
than the isotropic case. In the limit ^ 00 (no perpen¬ 

dicular pressure) we have X = (f’’ + \/( 7(8 -|- 9 ( 7))/(4 -I- dn). 
X —> !y/2 + a jA as (7 —>■ 0 and x 1 ~ 2/(30") as o" —>■ 00 . 
In this case, Jjx/d'p = 2 /( 72 A) ^ 1.27, and the minimum 
occurs at o" = 0, = 7 r/ 2 . Although parallel pressure domi¬ 

nance could alleviate the tension between the deflection and 
size, it is difficult to imagine a mechanism that would pro¬ 
duce mostly transverse random motion of the particles just 
behind a shock. 

The third possibility is that the upstream flow is a 
striped wind. In Append ix | A2| we generalize the shock jump 
condition of |Lyubarsky| ( |2003[ ) to include a possible 9 com¬ 
ponent of upstream magnetic field and an oblique ahead of 
shock velocity. We find that the shock jump condition is un¬ 
changed if we introduce an effective cr that takes into account 
the reconnection of alternating magnetic flux, in agreement 
with Lyubarsky (20031. 


3 CURVATURE PROBLEM 

In the simple shock model the emission is supposed to 
arise from a very short length behind the shock and can 
be approximated as a surface emission. This is illustrated 
in Figure (a), where the shock surface has a toroidal 
form. The emissivity at the shock surface is (see jv,n oc 


Q2+a^B' sin'c< 7 ')^+“, where & = 1/[7(1 — /3 • n)] = 1 /[7(1 — 
P cos Ip)] is the outflow Doppler factor, ip is the angle be¬ 
tween the line of sight and the flow velocity, B' and ro' are 
the magnetic field and the particle pitch angle in the fluid 
rest frame, respectively. The sinvu' factor only makes the 
feature more “boxy” (see (5.31 and the curvature is largely 
determined by the Doppler factor. Assuming axisymmetry, 
7 and B' are only functions of the polar angle 9, and ip 
is also dictated by the 9 variation of the flow velocity. (It 
would be isotropic if we were looking at a spherical surface 
with radial velocity everywhere.) It is then easily seen that 
the emissivity will roughly follow the projected constant 9 
contours, which are determined by the torus radius and its 
inclination. (A numerical illustration of the emissivity at a 
typical shock surface with artificial shape is shown in Figure 
i) Therefore, if the knot shows any curvature, it should be 


convex toward the pulsar (a “frown”). However, Rudy et al. 
| |2015 1 showed that in both Hubble and Keck images the arc 
corresponding to the knot is concave towards the pulsar (a 
“smile”). Thus, the observed curvature seems to contradict 
the simple shock model. 


4 DOPPLER-DEPOLARIZATION PROBLEM 


If we adopt the best-guess spectral index, a ~ 0.8 (Soller- 
|man|20d^ [Melatos et al.|2005[ ), the associated degree of lin¬ 
ear polarization for a uniform fields is ~ 0.73 (e.g., Rybicki 
fc Lightmanjl986 1 . This applies in the fluid rest frame. How¬ 
ever, we observe the sum of Doppler-boosted photons from 
different parts of the source. Now when a single photon is 
aberrated into our line of sight, its electric and magnetic 
vectors rotate in space about an axis along k x v and the 
observed electric vectors will not all be parallel and this con¬ 
tributes a degree of Doppler-depolarization. Typically the 
integrated degree of polarization will be ~ 0.4 — 0.5 (e.g., 
Lyutikov et al.|[2003 l. Special conditions are needed to re- 
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X 

(a) (b) 

Figure 3. An example showing the emissivity on the surface of the shock, (a) The 3D geometry of the shock surface. In this particular 
example, we assume the shock surface is given by r = sin^ 6 and the viewing angle is 9^^, = 2'n-l3, while the magnetization (7 = 1 and 
the spectral index a = 0.8. The contours correspond to Doppler factor ^ on the shock surface (the innermost contour has ^ = 8 and 
the interval between successive contours is 1). The location of the pulsar is shown by the red cross and the arrow indicates the pulsar 
rotation axis, (b) The contour plot of emissivity near the presumed site of the knot, as projected onto the plane of the sky, overlaid with 
constant 9 curves (black, dashed lines). The location of the pulsar is indicated by the black cross. We assumed that B oc Ijr and particle 
density n o(. Ijr^ ahead of the shock. We used the shock jump condition in Equations The procedure of calculating the emissivity 

is outlined in §5.2| 



Figure 4. Geometry of the shock front. For Crab Nebula the 
viewing angle should be 9ah = 27r/3, but here for better presen¬ 
tation we’ve plotted the case Squ < tt/2. The line of sight is along 
n and the flow velocity is denoted by v. Magnetic field is toroidal. 


cover the observed polarization of ~ 0.6. This will be dis¬ 
cussed in more detail in jj5.5| 


5 DETAILED SHOCK MODEL 

Having identified three problems with the simple shock sce¬ 
nario, we now turn to more detailed models to seek ways 
to address these problems. In this section we scrutinize the 
emission from the shock and derive the intensity map self- 
consistently, taking into account post shock flow dynamics. 
This gives us insights into limitations of the shock model. 


and suggests what other elements may be needed in a better 
model. 

The shock geometry is shown in Figure Since the 
knot is quite compact, we anticipate that the post shock 
flow Lorentz factor is large enough for us to expand about 
the center of the knot. 

We first set up a spherical coordinate system with the 
2 axis along the symmetry axis of the shock (the pulsar 
rotation axis), and write down the emissivity as a function 
of the spherical coordinates {r,0,(j>}. The direction of the 
line of sight, denoted by the unit vector n, is on the x — z 
plane, with Cartesian coordinates n = (sin0ob,0,cos0ob)- 
(We use {} for spherical coordinates and () for Cartesian 
coordinates.) For the Crab Nebula, we set Ooh = 27r/3 (e.g., 
Ng & Romani 2004 Weisskopf et al. 20121. The magnetic 


field B is assumed to lie along the (j) direction. 

The particle distribution in the fluid rest frame is sup¬ 
posed to be an isotropic, power-law distribution: 

f'{Y,zj')dYdQ.'dV' = (6) 

where zu' is the pitch angle of the particles in the fluid rest 
frame. If the outflow bulk Lorentz factor is 7 , the emissivity 
of the outflow would be 

>,n(i', n) = = i//^,n') = 


(7) 


where 


^ = 1/[7(1 - • n)] = 1 /[ 7(1 - ^ cos V>)] ( 8 ) 

is the outflow Doppler factor, tp is the angle between the line 
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Implications of Crab Inner Knot 5 


of sight and the flow velocity, and 


s a + 5/3 %/3 a l\_/a 5\ f 27rmc\ 

The intensity is then 

h,n{i', = J n)ds, 


omit the subscripts v and il on j or / in the following cal¬ 
culations but readers should keep in mind that these refer 
to quantities per steradian and per unit frequency.) 

5.1 Post shock flow dynamics and emissivity 

The flow should be axisymmetric and poloidal. If we consider 
a small bundle of flow lines that has a cross section 5A = 
51 ■ rScp, where S£ is the width measured in the poloidal 
plane and rSifi is that measured in the toroidal direction, for 
a steady flow we have the following conservation laws: 

• Conservation of particle number: nu ■ SA = const; 

• Conservation of energy: [{nm + + '^^]SA = 

const; 

• Conservation of magnetic flux (or V x E = 0): vb- S£ = 
const; 

• Equation of state: P = kK, where fc is a constant for 
this flow line. 

If the outflow were spherically symmetric and the post shock 
Lorentz factor is sufficiently high, B' oc 1/ur and n oc Ijur^, 
and if the particles evolve adiabatically, k' oc so 

2+a 

1 OC ' 


^-{5c,/3+2)^-(7c/3+3)^ 


7(1 — Pcostp) 

In the large 7 limit, we have 

j oc (1 +^2^2)-(2+“)^-(2+“)y-2-/3^-(W3+3) 


( 11 ) 


( 12 ) 


COS« 1 — ^ sin^ 6 'ob</>^. 


(14) 


72(1 - n • / 32 ) 


(15) 


Since the magnetic field in the fluid rest frame is along cj> di¬ 
rection, B' = (— sin (j), cos (f). Oh the angle between the pho- 
i 2 ~“ton and B' is (e.g., Lyutikov et al.|2003|: 


and the line integral in Equation (10 1 would be dominated 
by emissivity close to the shock surface. We consider the 
influence of curvature in the flow below. 


5.2 Emissivity at the shock front 

Suppose that the shock surface is a smooth function of 6 : r = 
r{6). The radial wind that goes through the point {r{9), 9, (f>} 
on the shock surface is deflected by an angle A and now is 
directed along a different polar angle 9-). Then t he angle tp 
between n and /32 is (e.g., Lyutikov et al.|2003 l 

cos Tp = cos 9ob cos 92 -f sin 9ob sin 92 cos p. (13) 


Since the radiation is mostly concentrated along the line of 
sight, i.e. 1^2 — fi'obi 1 , we can expand around the flow line 
that is directed toward us. To lowest order approximation, 
we have 


A photon propagating along the direction n in the nebula 
frame is emitted along the direction of n' in the fluid rest 
frame (e.g., Lyutikov et al.|[2003 l: 

, _ n-b 72 ^ 2 ( 7^11 1 ) 


( 9 ) 

COS vj' = YV • B' = 

so 

n - B' — sin 0 ob sin p 

(10) 

72(l-n-/32) 72(1 

-n-/ 32 )’ 


^ sin^ 0ob sin^ tp 

1/2 

(We 

sin tier = 

7|(1 —/J2 cosi/)^. 



(16) 


(17) 


In the following we derive 02 as a function of 0. At a 
point {r( 0 ), 0 , (f>} on the shock surface, the incident angle 5i 
is 


cotdi = 

r d9 


(18) 


and the curvature of the shock in the meridional plane is 


de^ 


According to the shock jump condition, the deflection angle 


(19) 


A = — arctan(y tandi). 


( 20 ) 


where y is the compression ratio. Consider the flow line 
that is directed toward us: 02 = 0 ob; writing this point 
as {r(0o), 00 , 0}, we have 0o -I- Ao = 0ob. Now if we devi¬ 
ate from this point by a small polar angle d9, we find that 
d 02 = d9 + dA, and since 

X sec^ 


dA_ 

dSi 1 + X^ tan^ di ’ 


dSi 

~d9 


/cir \2 _ , 

\deJ 


de2 




dA dA dSi 
d 0 ^ dfihB' 


( 21 ) 


( 22 ) 


(23) 


we get 


Ad I ^03 ^2 /.. COS 52 , , _ /o/ 1 \ 

d02 = d0 X-y^ + (l-X - =ad9, (24) 


where 


a = x- 


" cos^ (5i 
cos^ 52 


+ 1 -X 


cos 2 ' sindi _ 
COS^ 52 \ Tk 


(25) 


" cos^ (5i ' V ^ cos^ 5i J Rc sin 5i 
The same result can be derived based on a geometrical 
argument, see Figure Suppose that the shock radius of 
curvature in the meridional plane is Rc- Consider a small 
portion of the incoming wind with an opening angle d9, the 
difference in upstream incident angle is 

d5i =—d9 + d9Krk/sinSi. (26) 

Since the shock jump condition gives tan <52 = X ts-ii < 51 , we 
have 


d52 = d5ix cos^ 52 / cos^ 5i. 

Thus the outflow opening angle is 

d 02 = —d 52 + d9rk/(Rc sin 5i) 

/ cos^( 52 , cos^52^ rk 

= X—y^ + (l-X- 


'"cos^ di 


^cos^ 5i' Rc sin di 


(27) 


= a d9. 
(28) 
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Figure 5. A meridional cross section of the shock. 


Pulsar 



we have, to lowest order o^ X, Y, 

272 272 




I + 72V l + Y2Xyxl + a^Y2YVyr 


( 32 ) 


1 - 


1 - 


4sin^ 6 lob 72 <(> 


1/2 


(1 +7|7/)2)2 


(l+7|X2/a;2+a27|y2/y2)2j 


1/2 


(33) 


« C{u)k'B'^+’^ 


272 

2 +a ■ 

_1 + 7|i/)2 



1 - 


272 


4sin^ 6 obl 2 4 ‘'^ 
(1 + 7|i/'2)2 

-I 2 +a 


1 - 


[i + Y2xyxi + a^Y2YyyU 


^yix^lxl 


(l+7|X2/2;2+a2^2y2/j^2)2 


2 (l+<^) 


, (34) 


where xq = sin( 6 lob - Ao)/sin0ob, Vo = rfe(cosAo - 
cot(5i sin Ao). If Aq < 1 , xq « r^, yo ~ rk{l — cot5i tan A); 
further more if the upstream velocity is quasi-parallel to the 
shock surface, namely —>■ 0 , we would have yo —>■ x^k ~ 
XXo, but if the flow is quasi-perpendicular (di —>■ 7r/2), we 
have yo ~ Xk ~ xq. This already hints that we are more 
likely to get a knot that is elongated perpendicular to the 
symmetry axis when the upstream flow is quasi-parallel and 
a is small as we now discuss. 


Figure 6. Coordinate system for the projected intensity map on 
the plane of the sky. The origin is located at the peak intensity 
point of the knot, corresponding to the flow line behind the shock 
that is directed toward us. The pulsar is located somewhere on 
the positive Y axis. 


Now consider the projection. If we expand around the 
flow line that is directed toward us, namely the point 
{r{0o), 9o, 0} in spherical coordinates or (xq, 0, Zq) in Carte¬ 
sian coordinates, we have 


dx 

dy 

dz 


dv 

— sin ^ + r cos 9 
d9 


Vk sm9od(fi, 


dv 

— cos 9 — r sin 9 
d9 


d9. 


e^Bo 


d9. 


0=00 


(29) 


The coordinate system we use for the image on the plane 
of the sky is shown in Figure The projected coordinates 
{X, Y) corresponds to a rotation of the original pulsar-based 
X, y, z coordinates 


dX = dy, 

dY = — cos 9ohdx -\- sin 9ohdz. 


(30) 


Since 9o + Aq = 9oh, we finally get 
X = Xk sin(6lob - AQ)d(j>, 

Y = —rj;(cos Aq — cot di sinAo)d0, 


(31) 


near the point where the flow line is directed toward us. 

Using all the relations derived above, assuming 72 1, 

neglecting variation of 72 , B' and k' over the emitting region. 


5.3 Intensity map 

More generally, the emissivity in the outflow is (in the fol¬ 
lowing we omit the subscript “ 2 ” and just use 7 , /3, u to 
denote the outflow Lorentz factor, velocity and proper ve¬ 
locity, respectively, not necessarily at the shock front) 


j ^ C(i/)fc'B'^+“ 


r 27 

2 + a 


4sin^ 9ah+4>^' 

1 -i- ++ 



(1 -|- ++Y 


(35) 

Now supposing that the outflow roughly follows a spherical 
expansion, we have B' oc 1/ux and n oc Ijux'^, and if the 
particles evolve adiabatically, k' oc , so 




1 - 


4sin^ 9obX^(j>^ 


(1 + 


(l +«)/2 


(36) 


where the constants have been absorbed into the new pref¬ 
actor 

The outflow may bend in poloidal direction. We intro¬ 
duce a parameter r = r* d9/dx to characterize this bending, 
so that at a radius r, the direction of the outflow has a polar 
angle 


9 ^ 92 s + T{x/xk - 1), (37) 


where 92 s is the velocity polar angle at the shock front. For 
simplicity, suppose r is more or less the same over the emit¬ 
ting region. 

If the outflow is not curved, then from the geometry, we 
can see that the angle ■(/> between the line of sight and the 
velocity at radius x satisfies (at small angle limit, valid when 


(l + a )/2 
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Implications of Crab Inner Knot 7 


a > 0 , namely the deflection is monotonic in the emitting 
patch) 


~ % sin^ {^23 - d^hf 

r-^ \r-rk + rk/aj 

~ + f - — r) i^s-eof (38) 

\r-rfe + rfc/ay 

_ X^rl Y'^rl 

^0^'^ Voir - rk +rk/ay ’ 


where subscript s denotes quantities at the shock front. 
Based on this, we can easily see that if the flow bends ac¬ 
cording to Equation (371, then 


r -2 

(2 ' k • 2 ^ (2 I 

'Ip ^ ^sm 0oh(p + 


2 
2^2 


Tk 


XKj 

XqT'^ 


+ - 


r -rk+ Tk/a 
Yrk 


(Oa - Oo) -I- T{rfrk - 1) 


yoir-Vk +rk/a) 


+ T{rjrk - 1 ) 


(39) 


T > 0 means bending toward the pole; r < 0 means bending 
toward the equator; r = 0 takes us back to the straight flow. 

Parameterizing the distance along the line of sight as s, 
we get the intensity from the line integral 


/(i/,n) 

H f 

J r. 


J dsj{v,\ 


dr 


i(i^, n) 


C: 




1 - 


4sin^ 0ah'y'^4>^ 


(1 -I- 


(l+a)/2 


(40) 


where ip^ is determined by Equation (39l. In the last line 
r • n is taken to be 1 because ip ~ 0 ( 1 / 7 ). The upper limit 
on the integral can be set to 00 since the emissivity is only 
significant close to the shock. 

Let’s look at a few examples. In the following we assume 
that the flow does not slow down rapidly so the change of u, 
7 and P over a short distance can be neglected. The integral 
in Equation (401 can then be carried out numerically. 


5.3.1 Straight outflow 

For straight outflow, there are several interesting points to 
make: 

• If we ignore the pitch angle dependence, namely, omit 
the sin“ tu' factor in Equation ([^, from the Doppler factor 
we get the following: 

- The width of the knot perpendicular to the symmetry 
axis is rj± ~ 2 rj;/ 72 . For example, if a = 0 and Y = 0, 
the integral can be carried out to give I(v) = C 2 (i 2 )(l -I- 
72 ^ 3 / 3 : 0 )”^. The intensity profile depends on a\ it is more 
peaked for a steeper spectral index, as shown in Figure 
(a). 

- The ratio between the widths of the knot perpen¬ 
dicular and parallel to the symmetry axis y±/rj\\ can be 
calculated as axojyo = asin(0ob — Ao)/[sin 6 ^ob(cos Ao — 
cot5i sin Ao)] = h. Observationally we have 7 u/ri|| « 2. A 
contour plot of the knot intensity map is shown in Fig¬ 
ure]^ (a) for the case 6 = 2. The parameter 6 depends on 
the shock radius of curvature Rc, incident angle 61 and a. 
Figurej^shows 6 as a function of di for different a and Rc. 


For the case of quasi-parallel incident (61 —>■ 0), 6 ~ a/x 
and the requirement on Rc becomes the following if we 
want » 7 ±/? 7 || ~ 2 

R ^fc(l-x) . rk{l-x) 
(^“l)xsm5i xsindi 

The trend is that for very high cr flows, one needs very 
small 61 and/or Rc to get large enough b. This seems to 
be another difficulty associated with a very high a model. 

• When we take into account the pitch angle dependence 
of synchrotron power, r]±_ decreases because the synchrotron 
emissivity depends upon the perpendicular component of 
the magnetic field. This effect is shown in Figure [^(b) and 
Figure[^(b). For a spectral index a = 0.8 — 1, r]±_ can become 
as small as 0.4 times the value estimated in ^and this helps 
to reduce ryx/Tp. However, is not reduced by this effect; 
in order to still get the right aspect ratio /^n, we would 
require 6 to be 2-3 factors larger than that given in the last 
paragraph. 


5.3.2 Curved outflow 


For a curved outflow as introduced in Equation (371, a few 


examples are shown in Figure [T0| It can be seen that when 
the flow bends toward the equator, the knot is more ex¬ 
tended on the side away from the pulsar; in contrast, if the 
flow bends toward the pole, it is the side close to the pulsar 
that becomes more extended. The effect of pitch angle de¬ 
pendence is still to contract the emission in the x-direction. 

This effect may be important in explaining the curva¬ 
ture of the image shown in Keck and Hubble observations. 
If we imagine the outflow could bend through an angle 1 /72 
over a distance < r*,. Figure 10 (a) shows that bending to¬ 


ward the equator could have the potential to give the right 
curvature in the image. The effect can be tested using more 
realistic models of the post-shock flow. 


5.4 Anisotropic particle distribution 

Imagine that downstream of the shock the particle distri¬ 
bution is anisotropic and follows double adiabatic evolu¬ 
tion, then the intensity profile would be different from the 
isotropic case above. One complication here is that, as the 
plasma element moves along the flow line, the pitch angle 
distribution also evolves. For simplicity, assume that the par¬ 
ticles conserve their two adiabatic invariants so that the per¬ 
pendicular and parallel momentum satisfy px oc oc 

5^/3 and py oc L~^ oc nB~^. The first relation is eas¬ 
ily justified as the particle Larmor radius is much smaller 
than the length scale of the shock, while the second requires 
sufficiently low collision rate as particles move along the 
toroidal magnetic field. When this is the case, the parallel 
and perpendicular pressure satisfy d/dr(P||i?n“^) = 0 and 
d/dT{P±n~^ = 0. If as before B oc r~^ and n oc r~^, 
we would have P\\ oc r~^ and Px oc —Py decreases 

faster than P±. A proper treatment should evolve the dis¬ 
tribution function along the flow line. Here, for simplicity 
and illustration, we assume that the emission is dominated 
by a thin layer right at the shock so it is instructive to look 
at the emissivity at the shock front. 
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72^X0 


Figure 7. (a) Intensity profile along the line Y = 0 for different spectral indices o=0, 0.5, 1, when the outflow is straight and the pitch 
angle dependence is not taken into account. The profiles are normalized at the peak intensity, (b) Same as (a), except that the pitch 
angle dependence of synchrotron power is included. 



- 2-1012 - 2-1012 - 2-1012 

rzX/xo r2X/xo r2X/xo 

(a) (b) (c) 


Figure 8 . Knot shape and intensity profile for the case of straight outflow, in the coordinate system of Figure]^ (a) Contour Plot for a 
special example with aspect ratio parameter 6 = 2 and spectral index a = 0.8, when the pitch angle dependence is neglected; (b) Same as 
(a), except that the pitch angle dependence of synchrotron power is included: the knot gets squeezed in the x direction; (c) same as (b), 
except that 6 = 3. The intensity is normalized by the peak value and the contours go from 0.1 to 0.9 with increment 0.1. Same below. 



r2X/xo r2X/xo r2X/xo 

(a) (b) (c) 

Figure 10. Knot shape and intensity profile for the case of curved outflow, in the coordinate system of Figure]^ The pulsar should be 
located somewhere on the positive y axis, (a) Contour Plot for the case 6 = 3, q = 0.8, 'y 2 T = —1/2 (namely, the flow is bending toward 
the equator), when the pitch angle dependence is neglected, (b) Same as (a), except that 72 T = 1/2 (namely, the flow is bending toward 
the pole), (c) Same as (a), but with the pitch angle dependence of synchrotron power included. 
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(a) 


b 



b 



. Rc= 0.1 

. Rc=0.2 

- Rc=0.5 

- Rc=1 

- R„=1.5 

- Rc=2 


. Rc=0.1 

. Rc=0.2 

- Rc=0.5 

- Rc=1 

- R„=1.5 

- Rc=2 


. Rc=0.1 

. Rc=0.2 

- R„=0.5 

- Rc=1 

- R,=1.5 

- Rc=2 


Figure 9. The parameter b = axglyo as a function of the inci¬ 
dent angle <5i for different a and shock radius of curvature Rc ■ Rc 
is normalized to rj., namely the radius of the shock at the knot 
measured from the pulsar. Shaded regions correspond to physi¬ 
cally favored value of 6: 6 > 2. (a) tr = 5. (b) cr = 1. (c) cr = 0.1. 


Suppose the particle distribution function can be writ¬ 
ten as 

f' {-y' ,'cu')d')' d^l' dV' = k' g{s,in'Dj')dy' dO! dV', 

(42) 

where ^(sinro') describes the pitch angle distribution. We 
have the emissivity 

>.n(R,n) = 

= C{v)k'{B' sin g(sin zu'). 

When the downstream Lorentz factor 72 is large enough, the 
emissivity at the shock front can be approximated as 


j ^ C(u)k'B'^+°‘ 


272 

2 +q 

■ 4 sin^ 7 b 72 7 " 

[1 + 72VJ 


[ (1 + ylr)^ \ 


X 9 



4sin^ 6 ohyl(f>^'\ 

(1 + 72V)M ) 


(44) 


We show a few contour plots of the emissivity for different 
functional forms sr(sinrc 7 ') in Figure [TT| It can be easily seen 
that when P± dominates, the emission is more squeezed in 
the X direction and vice versa. (The parallel pressure domi¬ 
nant case may produce bipolar features.) If we focus on the 
former and assume g{sin-uj') = sin'^ti;', the width of the 
knot perpendicular to the symmetry axis can now be esti¬ 
mated as g± ~ 2rk -I- 3a -|- 25 according to Equation 

(441. In order to obtain « 0.5, we need (1) g ~ 15, 

namely, the pitch angle distribution is highly concentrated 
near 90°, (2) a ~ 0 and (3) di ~ 0. At the same time, to 
keep the right aspect ratio g±/g\\ ~ 2 , the radius of curva¬ 
ture of the shock surface should satisfy Rc < 0 . 2 rj;/sin di, 
while Tfe = 'kpD/A should be less than Teq. These strong 
requirements make the knot to be a very special point on 
the shock surface, and it’s doubtful whether such conditions 
can be always satisfied as the shock moves around due to the 
change of nebula pressure. The large perpendicular pressure 
is also not generic when there’s sufficient scattering in the 
flow. This remains hard to understand. 


5.5 Polarization 


In this subsection we examine in detail the Doppler- 
depolarization effect and suggest possible solutions. Let’s 
consider a plasma element with spherical coordinates 
{r{0),6,(j>} at the shock front. We follow the notations of 
Section |5.2| The emission has linear polarization perpen¬ 
dicular to the magnetic field in the comoving frame of the 
plasma, namely, e' — n' x B'. Upon a Lorentz transforma¬ 
tion back to the nebula frame, the direction of the electric 
held becomes (Lyutikov et al.|[2003l 


n X (B -I- n X (/)2 X B)) 

y'(l-n./3>-(B-11)771 


(45) 


since B ■ P 2 = 0 in our case. Introduce I = y which is a 
constant unit vector perpendicular to the plane containing 
n and z, the polarization position angle x is determined by 


cosx = e • (n X Z) 

_ cos 4 >{i — P 2 cos Ip) — P 2 sin 62 sin Oos si7 <p 
7(1 — 7 cos ip)'^ — sin^ Ooh sin^ P/li 
sin X = e ■ I 

sin (p(P 2 cos 62 — cos 7b) 

-7(1^^77a5s7P^-7r?7)b7in77^ 

Still assume that particle distribution satishes Equa¬ 
tion § and let p = 2a -I- 1 in the following. We consider 
steady state and suppose that radiation is mostly emitted by 
a thin shell with thickness Sr right after the shock. The lat¬ 
ter assumption is justihed if the how does not slow down too 
rapidly and roughly follows spherical expansion, as shown in 
§5.1| Now we can calculate the Stokes parameters by inte¬ 
grating over the visible surface of the shock front: 

j _ P + f d(p f d9 sinO r^/r'^ + (dr/dO)^ 

P + 1 D Jq Jq 

X (48) 


(46) 

(47) 
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- 2-1012 - 2-1012 - 2-1012 

yiX-lxo yiX-lxo 72^1X0 

(a) (b) (c) 


Figure 11. Contour plot of the emissivity on the shock surface, as projected on the plane of the sky, for different functional forms of 
particle pitch angle distribution, according to Equation ( |44| l. In these plots we take b = 3, a = 0.8. (a) Isotropic case, namely, g = 1. (b) 
^'(sintzj') = sin^ zu' with g* = 4. (c) g'(sin^') = (1 — sin^ with q = 1. 


r n27r ptt 

Q = A{u)j^ / dcj) d9 sm6 r\/r^ + {dr/dOY 
^ Jo Jo 

X (49) 

r p27t ptt 

U = A(u)-p— / dcj) d9 sin 9 r\/r^ + (dr/d9)^ 

D Jo Jo 

X (50) 

l/ = 0, (51) 


where 


A{u) 



3p—l\p/3p + 7\ e® /2TTmc\ 
12 ) V 12 ) mcP V 3e ) 


Since U integrates out to 0 under axisymmetry, the 
polarization degree is 11 = Q/I. We can get some analyt¬ 
ical results under ultrarelativistic limit. In this case due 
to the Doppler factor only a small region 0 ~ I/ 72 , 
02 — 9ob ~ 1/72 contributes significantly to the integral, so 
we can allow the integral limits to extend from —00 to -too. 
Also assume 72 itself as well as other quantities, Bi and 
niui, do not change significantly over this small patch. We 
change the variables to ^ = 72 ^ and rj = 72(^2 — ^ob), car¬ 
rying out the expansion in the limit I /72 <C 1, C ~ 0(1), 
rj ~ 0 ( 1 ), we have 




272 


1 - 


cos 2 x ■ 


1 + sin^ 9ob ’ 

4^^ sin^ 6 >ob 

(1-I--I-sin^ 6lob)^. 

(! + ??"- e" sin" e^b? - we sin" ^ob 

{1 + rj^ + sin^ 9ob)'^ - 4:^'^ sin^ 9ob 


(52) 

(53) 

(54) 


Also note that d92 = a d9 where a = x cos^ <52 / cos^ (5i + 
(1 — xcos^ < 52 /cos^ 5i)rfc/(sin5ii?c), so the integral over 9 
can be written into an integral over 6 ^ 2 . Let = ^sin 6 ^ob- 
When 9oh 1 / 72 , the result of polarization degree no longer 


depends on 9ob- 

p+1 


--Y- 


P + 7/3 


f r+oo r+00 

1/ 

I «/ —00 J —00 


X [(1 + rj^ + fif - 4e?]^’’“"’/^[(l + p^- - 47j\l] 


/ + 00 ^ + 00 

-00 J —00 


(p+l)/4 


(55) 


This gives the same result as Lyutikov et al. (20031. When 


p = 3, n = 9/16 = 56.25%, and when p = 2, 11 = 43.4%. For 
the oblique termination shock of the Crab Nebula, we expect 
the particle distribution power law index p to be between 2 

3, thus the theoretical upper limit for the polarization 
degree, under ultrarelativistic assumption, would be < 56%. 
In particular, for the best-guess spectral index a = 0.8, p — 
2.6, we have 11 = 51.6%. 

However, the flow may only be mildly relativistic, so the 
aberration would also be modest. In this case a more careful 
treatment is necessary and the result should depend on the 
details of the shock geometry. As an example we consider the 
hducial shock shape r = sin"^ 9 as shown in Figure where 
the downstream flow is only mildly relativistic, with Lorentz 
factor 72 ~ 3 — 5 at the knot (see Figure We evalu¬ 
ate the Stokes parameters I and Q in Equations (48 1 (491 


by integrating the emissivity j over the surface area where 
j ^ ii/, peak/10: we caution that this size is much larger than 
the observed value (see Figure]^), as expected. Neverthe¬ 
less, we proceed with the polarization calculation. We find 
that when p = 3, H = 68.9%; when p = 2, H = 61.2%; when 
p = 2.6, n = 66.2%. The results seem to be consistent with 
the observed polarization degree, but as we already noted, 
this is an artificial example that does not give the right com¬ 
pactness and curvature. To solve the polarization problem 
we need to understand the other two puzzles simultaneously. 

Clearly more careful modeling of the post-shock flow 
and its emission is needed. As the polarization degree is 
sensitive to the spectral index and the flow Lorentz factor, 
we could in principle constrain them from observed polar¬ 
ization degree. Observationally, it would be very helpful to 
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Implications of Crab Inner Knot 11 


get an accurate, simultaneous measurement of the spectral 
index. 


6 FLUXES OF THE KNOT IN DIFFERENT 
WAVE BANDS 

Having explored the challenges to the shock model in terms 
of the compactness, curvature and polarization of the knot, 
we now turn to other signature properties of the knot that 
can equally set important constraints on general shock mod¬ 
els. The first one among these is the flux of the knot in dif¬ 
ferent wavebands, which indicates the amount of dissipation 
in the relativistic outflow. It is commonly presumed that 
the termination shock is responsible for accelerating most 
of the IR- to X-ray-emitting particles, but it’s not yet well 
understood what mechanism could act efficiently enough at 
such a relativistic, perpendicular shock. Another very inter¬ 
esting question is whether the shock is the sole source for 
accelerating IR/optical emitting particles, or some in situ 
acceleration further out in the body of the nebula is needed. 
In this section we try to get some hints to these questions 
by characterizing the radiation efficiency of the knot and 
deducing the energy composition of the shock outflow. 


6.1 Observations 


Currently we have flux data in IR, optical and upper limits 
in radio, 7 -ray bands: 


• IR (Keck): Weisskopf et al. (20131 give a typical K' 
band knot magnitude K' = 15.9, after dereddening with the 
best fit value E{B — V ) = 0.52 ([Sollerman et al.||2000[ ) and 
the extinction law from |F itzpatrick (19991, this corresponds 


to FC°*' = 3.1 X 10“^^ 


ergs ^cm ^Hz 

,38 


. We define e{v) = 


vLiy^o,! (E j where E = 5 x 10^® ergs“^ is the pulsar spin 
down power and — D^EC°^ is the knot luminosity 

per steradian, as the radiation efficiency of the shock along 
the line of sight in the particular waveband, and so e(1.4 x 
lO^'^Hz) = 4.2 X 10"^ _ 


The latest Keck AO observations by Rudy et al.| (20151 
give the knot-pulsar flux ratio ~ 0.063 on average. If 
we take the photometry of the pulsar from |Sandbe^ 
& Sollerman (20091, the pulsar Ks band magnitude is 


13.80 ± 0.01, this corresponds to a dereddened flux of 
2.2 X 10“^®ergs“'^cm“^Hz“^, giving a knot flux in the K' 
band = 1.4 x 10"^Vgs"^cm-2Hz-^ 

• Optical (Hubble): The latest analyses of HST monitor¬ 
ing of the Crab (Rudy et al.|2015 I give the dereddened green 




a radiation efficiency of e(5.4 x 10^'* Hz) = 5 x 10 similar 
to the infrared estimate. 


• Radio (VLA ): The latest VLA observation (Bieten- 
holz et al. 20151 set an upper limit eC°^ < 300/rJy over 
a beam size of 1" x 1". Thus the radiation efficiency is 
e(5.5 X 10® Hz) < 1.6 x 10"“. 

• 7 -ray (Fermi): Although there isn’t enough spatial reso¬ 
lution to detect structures within the nebula, the total syn¬ 
chrotron 7 -ray flux of the whole nebula in quiescent state 
gives an upper limit on any possible synchrotron 7 -ray emis¬ 
sion from the knot. According to Buehler et al. (20121, 


the quiescent nebula synchrotron emission has an integral 


photon flux above 100 MeV of (6.1 ± 0.2) x 10“^cm~"s 
and a photon index of 3.59 ± 0.07. This corresponds to 
an energy flux above 100 MeV F{i> > 2.4 x 10®® Hz) = 
1.6 X 10“^® ergcm“® s“'^. So the upper limit on the radia¬ 
tion efficiency of the knot in 7 -rays is e{u > 2.4 x 10®® Hz) < 
1.2 X 10"®. 


6.2 Emissivity and particle enthalpy fraction 

For a broad spectrum of electrons/positrons, the pressure 
will be dominated by the part of the distribution that has 
a power law index of 2 , which corresponds to an emission 
spectrum F^, oc locally. A few spectral measurements 


have been done so far, for example, Sollerman (20031 and 
[Melatos et al.| ( |2005| ) give an index 0.8 from IR to NUV. 
Some more recent ones show slightly different spectral in¬ 
dex ranging from 0.63 (Tziamtzis et al.|2009l to 1.3 (Sand 


|berg &: Sollerm^|2009 1 but unfortunately these are based 


on non-simultaneous data. The spectrum of the nebula has 
an index ~ 0.6 — 0.8 in the same band. If the knot index 
is indeed of similar value, the IR/optical emitting particles 
should contribute a major fraction of particle pressure in the 
shock outflow. 

We estimate the particle enthalpy fraction in relevant 
energy bands from the radiation efficiency derived above. 
Suppose that the emitting area is A ~ (r'fe/ 72 )®, and the 
emissivity drops off significantly over a distance r*,/?, then 
we have 

~ vj,.,ciArkC = ( 4 ) C ~ 

(56) 

Assuming an isotropic particle distribution in the fluid rest 
frame, we obtain 




1 „ / 2 B'" 

-n2a.C7 ^ 


72.7 me 
2 t' 


(57) 


where ~ c ?/synchrotron cooling 

time of the particles in the fluid rest frame. Now the enthalpy 
fraction of the particles that emit at a particular photon 
frequency v is 


, gUy me -f 2 rkC 

vC ) = ' —• 

E/ (47r) 


(58) 


Using the above results we get 

/ ^ T2^cooi _ 8 72fcooi ^ 

3 rfc/c EI{A-k) 3 ru/c 

From this we can also estimate the particle injection rate 
per unit steradian 


(59) 


iV(7') = 


E 


A'^'me? 47r72 

Taking a typical downstream magnetic field B 2 = 5-3 mG, 
Th = T\i cm, Lorentz factor 72 ~ 5, and pulsar spin down 
luminosity E = 5 x 10®® erg s“®, we obtain the particle 
injection rate and enthalpy fraction in different energy range 
as follows: 

• IR-emitting particles: 

iViR = 4.6 X 10®® s“®sr“® cr//El 3 ( 72 / 5 )®, (61) 

7?iR = 0.1U(72/5)®BZ|/®r7; (62) 
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Optical-emitting particles: 

No = 1.5 X 10^® s-^sr-\riV51^(72/5)", (63) 

no = 0.075 ?(72/5)®BZ|/ViV; (64) 


• Radio-emitting particles: 

Nr < 4.6 X 10®® s“®sr“® crZjl Big (72/5)^ 
nn < 7.2 X 10~^?(72/5)®R:|/®riV. 


(65) 

( 66 ) 


• 7-ray emitting particles that give synchrotron photon 
at 100 MeV: 

N-, < 1.6 X 10®® s“®sr“® crZ-®513(72/5)®, (67) 

n-y < 0.005 <;(72/5)®5Z|/®rr7®. (68) 

We know from observations that the total optical to X-ray 
emission of the whole Crab Nebula requires a particle injec¬ 
tion rate of iQ®® ® pairs/s and the radi o emission re quires 
an injection rate of 10^® pairs/s (e.g., Arons 20121. The 
above results indicate that the oblique part of the termina¬ 
tion shock cannot provide enough radio-emitting particles. 
The number of IR/optical-emitting particles also seems in¬ 
sufficient. The IR emitting particle fraction 771R ~ 0.11 is 
indeed substantial, indicating that they might be around 
the peak in the post shock particle distribution. 

In order to investigate the acceleration mechanism, we 
can estimate how many orbits the particles complete to give 
the observed radiation power. Since the flow time scale is 
tflow = Tk/Rc), and the gyro frequency of the particles is 
Z2g = 72i^g = eBI{2'KY'm) in the nebula frame, we get the 
number of orbits to be (take IR band for example, and let 
u = (21410®"® Hz) 




eB 


?c 2 'Km 
6 X 10^ 


3e5 

4nmv 


3 ^ 


.72 


9c 


n IGn^reV 


(69) 


o3/2 -1/2 
^-^14 n?. 


Therefore, IR/optical emitting particles could in principle 
be accelerated by mechanisms that operate over many gy¬ 
ration time scales. However, this would be difficult for 7-ray 
emitting particles. 


7 VARIABILITY 

Another salient feature of the knot is its time variability. 
Both HST and Keck observations on timescales of months 
to years ( Rudy et al.|2015 1 found that variations of the knot 
size n± and 7711, surface brightness Rv) and flux F(y) cor¬ 
relate with the knot-pulsar separation ’5p. This is probably 
not that surprising as the inner nebula is known for its vari¬ 
able emission, especially the dynamic wisps, ever since the 
pioneering work of e.g. |Scargle| ( |1969| ). Here we show that 
based on quite general grounds shock models can give a nat¬ 
ural explanation to the variability of the knot. 

In the shock model, variability can arise from at least 
two sources: changes in the surrounding nebula that result 
in modified boundary condition for the downstream flow, 
and variations in the upstream condition due to the pulsar. 


7.1 Variations initiated downstream 

Inhomogeneous and unsteady pressure in the nebula can 
cause the shock radius and shape to vary over time. Nu¬ 
merical MHD simulations (e.g. 2D runs by [Camus et al.| 
2009 and 3D runs by Forth et al. 20141 show that even if 


the upstream condition is fixed, the post-shock flow is highly 
variable as observed, due to vortex shedding from the ter¬ 
mination shock, and interaction with large amplitude waves 
in the nebula. The simulations of these complex variations 
exhibit timescales consistent with the observations. 

Suppose that some change in the nebula pressure causes 
the shock radius Vk and incident angle hi to change at the 
location of the knot. Assume the time scale to be long 
enough so that the shock evolves quasi-statically. Accord¬ 
ing to Equation Q, j(i/) oc k'B'^^°‘. Suppose k' oc 712, 
and since 722 ~ n\Ui/u 2 ~ nyuxjni (if 72 is large enough), 
B' — B 21^2 ~ 5 i/(x 72), at peak intensity point & ~ 272, 
we get j oc niU\{B\/xY'^°‘■ Also noticing that niwi oc r^®. 


Bi oc rI, , we find j oc 


-(3+a) 


Now suppose the emission 


length is £ (which may depend on too), then the peak in¬ 
tensity IR) oc Tj, The knot-pulsar separation can be 

written as = rk^/D, so it’s not surprising that IR) is 
correlated with tkp. For example, if £ oc rt and if A weren’t 
changed significantly, we would have IR) oc Also 

the width n± ~ i'k/"l 2 , is again correlated with and thus 
and roughly 77X oc Tp. Note that in addition to r^, these 
properties also depend on 72, which is essentially changing 
with as can be seen from Equation The variation of 
the width parallel to the symmetry axis n\\ is more involved 
as one needs to take into account the variation in the shock 
shape. [Forth et al.j ( [2014[ ) reported a correlation between 
knot flux and its separation from the pulsar that is similar 
to the observed one based on their 3D MHD simulation, sug¬ 
gesting that the shock model could give a good explanation 
of the knot variation. 


7.2 Variations initiated upstream 

Alternatively, a variation in the pulsar wind could cause 
variable knot emission. In such a scenario, the shortest pos¬ 
sible variation time scale is ~ — cos(l/72))/c ~ 

rfe/(27|c) ~ 1 day(rfc/10®®’ cm)(72/5)“®. Suppose the down¬ 
stream pressure is fixed, while the upstream energy flux 
in the direction along the line of sight changes. In such a 
case, the shock radius and shape would still vary accord¬ 
ingly, but as jR) oc fc'®®'''“5'®+“ oc 712 ^ 2 B^'^, it won’t 
change much if 72272 and B 2 stay more or less the same. 
However, the knot-pulsar separation 4tp = rkl^/D, and knot 
size 77X ~ rfe/72 would depend on both rk and 5i; as a result, 
while the correlation between n± and 'I'p is similar to that 


in f 7.1 the correlation between IR) and tkp would be quite 


different. 


8 DISCUSSIONS 


In this paper, we have presented a discussion of the impli¬ 
cations of recent observations of the inner knot of the Crab 
Nebula for its interpretation as a flow immediately behind 


a pulsar wind termination shock | Komissarov & Lyubarsky 
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2004 Komissarov fc Lyutikov|2011 1. Some of these observa¬ 
tions can be characterized as being broadly consistent with 
this interpretation; others present challenges to it. 

In the former category, we remark that some sort of 
feature is expected. Indeed, it is quite surprising that the 
radiative efficiency of the outflow of a supposedly ultrarela- 
tivistic wind in the spectral band where most of the pulsar 
power is radiated, at near UV frequency, is apparently so 
low (e < 10“®). (In principle, the wind could be more radia¬ 
tive in 7-rays, but the absence of a correlation of the knot 
variation with the 7-ray flares can be used to argue against 
this.) Another success of this model is that the knot is lo¬ 
cated on the projected nebular symmetry axis (presumably 
coincident with the pulsar spin axis), that it is elongated 
perpendicular to this axis and that the direction of the lin¬ 
ear polarization, supposedly perpendicular to the expected 
toroidal magnetic field, is along the symmetry axis. In prin¬ 
ciple, these attributes could be associated with a mildly rel¬ 
ativistic internal shock, or a Mach disk, in an inward ex¬ 
tension of the southern jet, but it is then surprising that 
the observed dissipation is seen so close to the pulsar and 
that the possible other knot observed with this position an¬ 
gle relative to the pulsar is too far away (~ 3.8”, see 


Hester 


et al. 19951. The variability on an outflow timescale is at 


least expected given the widespread, observed activity of the 
“wisps” in the inner nebula ( Scargle|1969 Hester et al. 1995 


2002). Furthermore the tentative correlation of the flux and 


location of the knot has a natural interpretation in terms of 
motion of the shock surface. 

The first challenge to this interpretation is that the 
compactness of the knot compared to its offset from the 
pulsar (?7 u/4'p ~ 0.5) and its elongation perpendicular to 
the symmetry axis ~ 2) were hard to reproduce 

quantitatively using the shock model. This may require the 
shock to have an almost grazing upstream velocity yet very 
small radius of curvature on the meridional crosssection 
[Re < 0.2rfc/sin5i, see (5.41, which seems unnatural, or re¬ 


quire the particle velocity distribution in the nebula frame 
to be elongated transversely, either through shear in the 
bulk velocity or a strong parallel pressure which is the op¬ 
posite of what would be expected at a perpendicular shock 
front. Even more perplexing is the curvature of the knot. 
It is concave towards, as opposed to away from, the pul¬ 
sar. As we have described, it may be possible to account for 
this but only with a special flow field or observer orienta¬ 
tion. Another great challenge is the degree of polarization. 
As we have explained even with a uniform toroidal field, 
the Doppler rotation effect will effectively reduce the syn¬ 
chrotron polarization of 0.73 to ~ 0.5, significantly less than 
the reported value. 

However, the greatest challenge to this model may come 
from understanding how the particle acceleration actually 
takes place. In models of relativistic cosmic plasmas, it is 
common to presume that particle acceleration takes place 
at relativistic shocks in much the same way as it does at 
non-relativistic shocks with particles gaining energy as they 
are scattered back and forth across the shock front. However, 
the kinematics of relativistic shocks are fundamentally dif¬ 
ferent and although transmitted distribution functions can 


be computed for given assumptions (e.g., [Gallant et al. 1992 


Hoshino et al.|lS)M Sironi &: Spitkovsky|2009 20111, a fully 

consistent description of the acceleration eludes us. There 


are two key questions that are raised by this interpretation. 
The first is why is the knot so strongly polarized if there are 
strong magnetic fluctuations close to the shock front and 
(j < 1. Somehow these fluctuations have to damp out leav¬ 
ing behind a dynamically weak yet very uniform toroidal 
held. The second question is whether or not the shock is the 
dominant source of particles with energies above, say, 100 
GeV emitting infrared synchrotron radiation or are most of 
the electrons and positrons accelerated much farther out in 
the nebula where most of the emission occurs. It is not pos¬ 
sible to give a dehnitive answer to this question but we can 


estimate how fast the ultra-violet- (Hennessy et al. 19921 
and X-ray- ( Weisskopf et al.|2000 Seward et al.|2006 1 emit¬ 
ting electrons would have to transport. If we suppose that 
the transport is essentially radially diffusive, then in order 
to reach the outer parts of the observed nebulae at these en¬ 
ergies, the mean-free path perpendicular to the presumed 
toroidal field direction, must be ~ 10^(73/10“"^ G)® gyro 
radii (independent of emission frequency). This is sufficiently 
large and may suggest that some of the particle acceleration 
is actually in situ. 

One of the longest debates involving the Crab Nebula 
is the value of the magnetization parameter, a. There is 
little doubt that the torque applied to the neutron star is 
essentially electromagnetic and the outward energy and an¬ 
gular momentum flux in the vicinity of the light cylinder 
is mostly electromagnetic; otherwise it would be very diffi¬ 


cult to account for the observed pulse proEles (e.g., Arons 


20091. However the wind flows through nine decades of ra¬ 


dius before its momentum flux matches the ambient nebula 
pressure and much can happen over this journey. Its initial 
electromagnetic structure is a combination of the axisym- 
metric, increasingly toroidal field associated with an aligned 
rotator and the striped, low latitude field reversing every 
~ 5 X 10® cm of an oblique rotator. The wind radial ve¬ 
locity is expected to accelerate to ultrarelativistic velocity 
- the velocity of the local Lorentz frame in which the local 
electric field vanishes. This can happen with minimal dissi¬ 
pation when the acceleration is due to the electromagnetic 
Lorentz force. However, non-stationarity at the source and 
instability in the outflow will inevitably heat the gas - we 
presume a pair plasma - and limit the acceleration of the 
outflow. This is especially important at low latitude where 
the stripes may be able to reconnect with considerable dissi¬ 
pation (e.g., Goroniti^l990 1 although strong arguments have 
been given that there is insufficient time in the comoving 
frame for this to be effective (e.g., Lyubarsky fc Kirk|2001 1. 

Consideration of the large scale current flow in the wind 
allows a complementary description to one based upon mag¬ 
netic field. On quite general grounds, the unipolar induction 
mechanism will generate a potential difference ~ 50 PV 
between the poles and the equator and there will be a 
quadrupolar poloidal current distribution, with a magni¬ 
tude ~500 TA in each hemisphere. The latitude-dependence 
of this current is what ultimately determines the shape of 
the shock surface. If the current “closes” well within the 
shock, the magnetization will be low and vice versa. When 
(7 is low, the ratio of the bulk acceleration to the “Ohmic” 
dissipation dictates the Mach number and asymptotic fluid 
Lorentz factor. Unfortunately we do not have a good enough 
understanding of these subtle issues to predict a and 7 and 
must rely on observational guidance. 
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We also have some clues from the nebula beyond the 
shock. What we observe is that electrons and positrons are 
accelerated to TeV energies and these particles ultimately 
radiate most of the rotational energy extracted from the neu¬ 
tron star. They do so on time scales ~ 30 yr., long compared 
with the flow time from the pulsar but slow compared with 
the lifetime of the nebula with high average efficiency al¬ 
though we know that energy can build up in the nebula and 


then be released faster than it is supplied (Wilson-Hodge 


et al. 20111. There is, of course, a broad electron distribu¬ 
tion from radio-emitting GeV electrons to the ~ 3 PeV 
7-ray synchrotron-radiating electrons that produce the dra¬ 
matic flares. (These also radiate ~ 0.01 of their power as 
inverse Compton 7-rays.) The low energy electrons accumu¬ 
lated over the nebula lifetime while the high energy electrons 
have to be accelerated locally as their cooling times are so 
short. At the highest energies they cool in a small fraction 
of a Larmor orbit. On these grounds, it has been argued 
that some of the nebula must be magnetically-dominated as 
this seems to be necessary to account for such dramatic, 
impulsive acceleration. However, it has also been argued 
( Atoyan &: Aharonian|[T9^ Meyer et al.||2010 l, on the ba¬ 
sis of the ratio of the Compton to the synchrotron power, 
that, on average, the particle energy density, dominated by 
TeV electrons, exceeds the magnetic energy density. This 
suggests that the nebula is quite inhomogeneous. The rel¬ 
atively high average nebula polarization in the optical and 
X-ray bands (~ 0.2, see, e.g. Woltjer|1957 [Port fc Walraven 


1956 Schmidt et al.|1979 Weisskopf et al.|1978| ) indicates a 

long range order again suggesting that f5 = SnP/B^ is not 
much greater than unity. 

Returning to the physical conditions at the inner knot, 
interpreted as a shock, we have argued that the compact¬ 
ness combined with the shape implies a low magnetization, 
a < 1. One way to reconcile this with the indication that 
some of the nebula is highly magnetized is if the line of sight 
passes through the striped wind at small radius and the dis¬ 


sipation of this magnetic energy is efficient (e.g. Komissarov 
2013 Forth et al. 20131. However, this seems inconsistent 


with the high polarization. We therefore consider some alter¬ 
native possibilities. The first of these is that the post shock 
flow is more complicated and typically contains a point of 
inflexion after it is pulled towards the poles by the magnetic 
tension. The path length along the flow will be long at these 
inflexions and the single flow line when the tangent is also 
directed towards us is the location of the knot. This can 
be examined using numerical simulations of the post-shock 
flow. The second idea is that the pre-shock flow is non-radial 
due to the action of magnetic stress and the knot is formed 
where a special path perhaps associated with a current sheet 
in the magnetosphere becomes tangent to the line of sight. 

The third possibility is more radical, that the outflow 
is highly magnetized and there is no shock. When most of 
the pressure and current is associated with the highest en¬ 
ergy particles that are radiating efficiently, a simple fluid 
description is called into question. The flow is better de¬ 
scribed using individual particle orbits and the deceleration 
is gradual. Well beyond the light cylinder the magnetic field 
is essentially toroidal and the poloidal current is composed of 
electrons and positrons undergoing gradient, curvature and 
magnetization drifts. It is possible that the details of these 
motions are affected by radiation reaction operating on the 


highest energy particles. As the guiding centers of the orbits 
move along the electric field, the particles gain energy and 
a self-consistent description of the particle acceleration and 
emission becomes possible. The nature of the inner knot is 
then associated with the non-radial motion of the guiding 
centers that leads to a maximum intensity on the sky. This 
basic idea can be explored using PIC codes. 

We note that [Lyutikov et al.| ( |2015[ ) have carried out 
a similar investigation of the knot emission as coming from 
the termination shock, and reached the conclusion that an 
effectively low a shock can produce all the observed prop¬ 
erties of the knot. They use the Kompaneets approximation 
where it is presumed that the post-shock flow is isobaric 
to model the shape of the shock; this only explores a quite 
limited parameter space as compared to our more general 
analytical expansion around the peak, while the peak in¬ 
tensity may be close to the boundary of the shock so that 
adding up the emissivity of the two folded sections of the 
shock may help produce a small enough 7711. However, the 
size perpendicular to the symmetry axis, rj±, still appears 
larger than the observed value in their model, and the cur¬ 
vature of the knot is inconsistent with current observation. 
[Lyutikov et al.| ( |2015| ) argue that hnite thickness of emitting 
region and/or limited resolution of the observation could 
change the curvature. They also point out that [Moran et ah] 
| |2013| ) used a very small aperture centered on the peak of 
the knot to measure the polarization; within such a small 
region the Doppler depolarization effect is not significant, 
especially when the flow is only mildly relativistic. But the 
high polarization can still be challenging in terms of parti¬ 
cle acceleration models that require turbulence at the shock. 
Undoubtedly, the inner knot of the Crab Nebula will be ob¬ 
served more and the conclusions on which we have based 
our analysis will be corroborated or modified. 

Despite the somewhat inconclusive nature of our inves¬ 
tigation, it is apparent that observations of AGN jets, 7-ray 
bursts and pulsar wind nebulae, especially the Crab Nebula, 
is rapidly enhancing our appreciation and understanding of 
how relativistic plasma behave in cosmic sources. Interest¬ 
ingly, similar knot structures are also seen in Chandra im¬ 
ages of Vela (Levenfish et al. 20131, the only other pulsar 
wind nebula which we see at comparable resolution. This 
may indicate the generality of some features and could give 
further hints. We believe future multiwavelength monitoring 
campaign of the Crab and other related sources will help us 
understand many of the aforementioned open questions. 
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APPENDIX A: GENERALIZED SHOCK JUMP 
CONDITION 

A1 Anisotropic plasma pressure 

In the general case when the particle distribution is gy- 
rotropic but anisotropic between the directions perpen¬ 
dicular and parallel to local magnetic field, the energy- 
stress tensor for the fluid part can be written as = 

[p -I- Pl) + P±g^'' + (P|| — P±) , where is a unit 

4-vector along the direction of magnetic field in the fluid rest 
frame. The energy density is p = nmc? -|- 2P± + Py. First 
consider a perpendicular shock where the incoming velocity 
is along the shock normal and magnetic field is perpendicu¬ 
lar to the shock normal. Applying conservation laws at the 
shock front gives 

niwi = n2U2, (Al) 

viBi = V 2 B 2 , (A2) 

1 0 1 9 

(pi +Plx)7lMl + = (P2 +P2x)72M2 + —VlBf, 

47r 47r 

(A3) 

(pi + Pix) Wi + Pix + (yi -f 1) Bi = 

(p2 + P2x) ih + P 2 X + ^ (nf + 1 ) Pf- (A4) 

Suppose that the upstream flow is cold (Pix = Pi|| = 0) 
and ultrarelativistic (ui « 71, iii ~ 1), and define a = 
Pi/(47rni7imc^ui). Let P211 = C-Pzx, then p2 = n 2 mY + 
(2 -I- C)P2X ~ (2 -f C)P2X under strong shock assumptions. 
Now the above system of equations can be solved to get the 
compression ratio 

X~V2 

2 -b 3(7 -b Ccr + V(2 + 5(7 -b 3C(7)2 -b 8(7(1 -b C)"' ("^5) 

■4(2 + 0(1 +(7) 

When = 1, we get back to the familiar result of isotropic 
case, Equation ([T|). In the limit C — 0 (no parallel pressure), 
we get X = (2 + 3(7 + \/4 + 28(7 + 25a'^)/{8 + 8(7). x 
1/2 + 3(7/4 when cr —> 0 and x 1 ~ 2/(5(7) when a —>■ 00. 
In another limit 00 (no perpendicular pressure) we have 
X = ((7 + 0'^(8 + 9 ct))/(4 + 4ct). X — 0n/\^+(7/4 as (7 —>• 0 
and X 1 “ 2/(3(7) as a ^ 00. Figure [AT] shows how the 
compression ratio x changes with for different a. If the 
upstream plasma is warm and anisotropic, we just need to 
change the definition of a into a = Pi/(47r(pi + Pix)7i)5 
where pi + Pix = u>i is the upstream enthalpy density. 

When the incoming velocity is oblique with respect to 
the shock front, it is easily shown that the component of 
velocity parallel to the shock surface is conserved so one can 
make a Lorentz transformation to get back to the perpendic¬ 
ular case. The relations for deflection angle and downstream 
Lorentz factor, Equations and Q, remain valid. 



Figure Al. Shock compression ratio x as a function of the 
anisotropy parameter ( = Py/Py, for different cr’s. 



Figure A2. Local coordinate system for deriving the shock jump 
condition in a striped wind. 


has both a 0 component and a (/> component, while its r 
component is zero. Consider isotropic plasma first: ~ 

Pg'^'' + (P + p)u^u'', and let the adiabatic index to be 
F = 4/3. For the upstream stress tensor, we average over 
a stripe wavelength. Continuity of the stress tensor compo¬ 
nents , field components B^, Ey, Ez and particle flux 
across the shock front give the following equations: 


A2 The upstream magnetic field has both a (f) and 
a 9 component, and is possibly striped 

We use a local coordinate system at the shock front, as 
shown in Figure |A2| where x corresponds to the shock 
normal, z corresponds to (j) direction and vi is in radial 
direction. We suppose that the magnetic field upstream 


pOx j B\y {^B\yV\x BlxVXy') + B^^'^lx 

477 


+ 71 {Pi + Pi) “la 


B2y {B2yV2x — B2xV2y) — B2z {B2xV2z — B2zV2x) 


47r 


+ 72 (P’2 + P2) U 2 x, 


(A6) 
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{BizViy + Bif) + — {BiyVix — BixViy)^^ 

+ (^-^{BizVix + Bi^viy + B\) + (Pi + pi) vfx + Pi^ 

= (P2 + P2) U2x + P2 

+ ^ [—2 {fB2zV2y — B2yV2z)^ + B^f) + {B2yV2x — B2xV2y)^ 

+ {B 2 x'V 2 z — B 2 zV 2 x)^ + (B 2 zV 2 y — B 2 yV 2 zf + Pf] ) 

(A7) 


: ^(Pl + pi) - 


B\xB\y 


471 


{B2xV2z — B2zV2x) {B2zV2y — B2yV2z) + B2xB2y 


471 


+ (P 2 + P 2 ) U2xU2yi 


(A8) 


/ l^lyPlz (Pli/77l3; Pla;l^ly) P PlrcPlz \ 

\ 4 ^ / 

— {B2yV2x — B2xV2y) {B2zV2y — B2yV2z) + B2xB2z 
471 

+ (P 2 + P 2 ) U2xU2z, 

(A9) 


{Bix) — B 2 x, (AlO) 

{,B\yi — p2yi 01 (Plzl^lx) — B2z^2x B2x^2z ^ (-^41) 

{Exz') — p2z) or IBxyVXx PlxUly) — B2yV2x B2xV2y ^ 

(A12) 

niWlz: = n2M2z: (^13) 


Here {) denotes an average over a wavelength. Suppose 
that the upstream flow is ultra-relativistic, namely ui « 1, 
7i « Ml S> 1. We can write Miz, = sint^i, viy = cos5i. From 
the geometry as shown in Figure we have Biz = B^, 
Bix = —Bg cos 5i, Biy — Be sin Si. The solution to the 
above equations satisfies V 2 z = 0, and the compression ratio 
is 


X = 


V2x 

'^Ix 


+ 0-2 + 14cr<; + 2 ct + 1 + ac + cr + 1 


6(ct + 1) 


where a is defined as 


_ {Bl) + {Bl) 


4717 ^ (Pi + pi) 


(A14) 


(A15) 


and c = ({Bgf + (B^) ") / {{Bj) + (P^)). if ^e let a = 
cr?/(CT + 1), we get 


X = ^ (l + cr + V+ 145- + , 


(A16) 


generalizing the results obtained by Lyubarsky (20031. 
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